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We analyze the impact of multichannel scattering in harmonic waveguides on the positions and widths of 
confinement-induced resonances for both isotropic and anisotropic transversal confinement. Multichannel scat- 
tering amplitudes and transmission coefficients are calculated and used to characterize the resonant behavior 
of atomic collisions with varying anisotropy. A mechanism is established which leads to a splitting of the 
confinement-induced resonance in the presence of anisotropy. 
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I. INTRODUCTION 



Confinement-induced resonances (CIRs), originally pre- 
dicted in the seminal work of Olshanii \^ and more recently 
observed experimentally for both bosons (0-31 and fermions 
151], have attracted a great deal of attention during the past few 
years. The immediate reason is that they represent a valuable 
tool for the control of the atomic interactions thereby allowing 
to enter and probe the regime of strongly correlated bosonic 
or fermionic many -body systems. Beyond this, changes of the 
transversal confinement potential possess direct consequences 
for the scattering behavior of the atoms in the waveguide such 
that a multitude of binary resonance profiles are accessible. 
Examples for the variety of scattering properties in waveg- 
uides are multi-channel confinement-induced resonances ja], 
resonant molecule-formation processes and center of mass 
coupling effects lul HD as well as a dual CIR l9J which leads 
to complete quantum suppression of scattering. 

A recent experiment lUOll on CIRs in transversally 
anisotropic waveguides has shown the necessity of an ade- 
quate theoretical approach for describing collisional processes 
in confined geometries of one- and two-dimensional char- 
acter. A main result of this experimental work lllOll was 
the observation of a splitting of the CIR-related loss sig- 
nals in the presence of a transversal anisotropy. Two re- 
cent attempts to explain this splitting have not been success- 
ful lllll[l2ll . The corresponding studies represent investiga- 
tions of the pseudojjotential scattering with a single open 
transverse channel 111] in the case of tightly confining waveg- 
uides. Both the width of the CIR and the full multichannel 
character of the problem were not taken into account. On 
the other hand it has been shown that the non-separability 
of the center of mass and internal motion for atomic scatter- 
ing in anharmonic (transversally isotropic) waveguide poten- 
tials changes the resonance picture qualitatively and new res- 
onances (anharmonicity-induced resonances (AIRs)) are oc- 
cur where molecular excited center of mass states cross the 
threshold ^ [T3[]. Very recently the coupling of the center 
of mass excitations in anharmonic isotropic and anisotropic 
confining potentials to the ground state has been analyzed in 
great detail in ab initio calculations 1 1411 . There a very good 
agreement of the AIR splitting with the distance between the 
maxima of the atomic loss in the experiment llOtl was found. 



However, it has not been analyzed what happens with the "har- 
monic" part of the investigated two-body spectrum in the trap 
and whether the "harmonic" CIR agrees with previous results 

In previous works ^ [III llJ, lo, llal t he specific ratio 
ai_/as ~ 1.4603... (where a^ = y^h/{jj.co±),as,IJ. and co± are 
the harmonic oscillator length, scattering length in free space, 
reduced atomic mass and harmonic oscillator frequency, re- 
spectively) yielding the position of the CIR was defined by 
the zero Im{f(){a±/as)} = of the scattering amplitude /o 
in the transversal ground state in the zero-energy limit. This 
point also corresponds to the absolute minimum of the trans- 
mission coefficient T(){a^/as) = |1 +/o(fl±/ai)P — >■ which 
represents an important scattering observable near the CIR 
Illl [lll[l2l[l5l[l6ll . In the present work, we explore a situation 
where not only the ground transversal channel but also excited 
channels contribute to the scattering process in the waveguide 
and demonstrate the importance of multichannel scattering ef- 
fects in the resonant region. 

Within our description of few-channel ultracold scattering 
in confined geometries we observe a mechanism which leads 
to the splitting of the CIR under the action of a transversal 
anisotropy of the waveguide. It is based on the fact that the 
total transmission coefficient T, which we define as a sum of 
partial coefficients 7] emerging from the different transverse 
(ground and excited) states labeled by /, 



Y^WiTi 



(1) 



averaged over the initial populations WJ, possesses its main 
contribution near the CIR from the first excited transversal 
state and not from the ground state. This follows from the ob- 
servation that To possesses a much deeper as well as broader 
transmission well around its minimum compared to the wells 
(minima) of Ti{i ^ 0) for the excited states near the CIR ]|a]. 
Thus, even for low-energy pair collisions the scattering prop- 
erties near the CIR location are not determined by the partial 
coefficient To, where Tq -^ 0, but by the behavior of the coef- 
ficients T;^ in excited states / 7^ 0. Employing this mech- 
anism we find a splitting of the minimum of the total trans- 
mission coefficient (1) in an anisotropic trap. The splitting 
is a consequence of the different dependencies of the partial 
coefficients Ti{a±/as) on a±/as for different states in the reso- 
nant region. A necessary prerequisite of the appearance of this 



splitting effect is therefore the occupation of excited transver- 
sal states in the waveguide. We emphasize, that the present 
investigation is performed for harmonic traps, i.e. we inves- 
tigate the influence of the anisotropy of a harmonic transver- 
sal confinement on the CIRs Ol [Til [H, [H, [l^ opposite to 
the above-mentioned works that explore the effects of anhar- 
monicity ||8l[T3|] and anisotropy lll4ll . 

In the following section we describe the computational ap- 
proach to our multichannel scattering problem in the confined 
quasi- ID geometry. The corresponding results are discussed 
in the third section. The short fourth section is devoted to a 
discussion of the mechanism leading in the harmonic waveg- 
uide with transverse anisotropy to the splitting of the CIRs. 
Finally we provide a brief conclusion. 




II. MULTICHANNEL SCATTERING PROBLEM IN 
ANISOTROPIC HARMONIC WAVEGUIDES 

To calculate the partial transmission coefficients (the single 
index / is here replaced by the double index ni,n2 indicating 
the quantum numbers belonging to the different transversal 
degrees of freedom of the waveguide) 
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describing the transmission probability from the initial trans- 
verse state («i,«2) to all possible final states («j,«2) in the 
course of the collision of identical bosons in a harmonic 
waveguide with the transverse trapping potential ^/^((B^x" + 
ofy^), we solve the multichannel scattering problem for the 
3D Hamiltonian 



H{x,y,z) = - — ^r + ^mU' + :^my' + V{r) , (3) 

depending on the relative variables r = (x,^,^), with the 
asymptotic scattering wave function for | z | ^ +°° 



V^«l,«2W =COs(/:„i.„2Z)0„i,n2(-^'3')+ Y^ fnl'Mi 
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FIG. 1. (color online) Dependence of the partial transmission coeffi- 
cients r„|.„, (axA'.s,£||) on a^/os and £11 in an isotropic waveguide 
0)1 /ob = 1. 



E^"^'"-' -^ e" = T'i(0±{n +1). The asymptotic wave function 
(4) is expUcitly symmetric with respect to the exchange of the 
atoms. 

In order to solve the multichannel scattering problem (3,4) 
we extend the approach developed in iflTl UM for scatter- 
ing in three dimensions: the expansion over the spherical 
harmonics on a grid (two-dimensional discrete-variable rep- 
resentation) is replaced by the expansion over the product 
states 0„j ,„(x,3') = (p„| {x)(p„2{y) where (p,,. are the eigenfunc- 
tions of the ID harmonic oscillator. The calculations are per- 
formed with the finite-range Gaussian approximation V{r) = 
— Voexpl— r^/rQJ for the interparticle interaction with Vo cho- 
sen such that we have one weakly-bound state in the potential 
well. We do expect that all effects shown below are to a large 
extent independent of the chosen interaction potential V{r). 
As an independent check, we have verified that the results 
obtained for the four open channel scattering in an isotropic 
waveguide Oi = 02 = C0± using the screened Coulomb po- 
tential la] are reproduced by the present approach. Note that 
for convergence typically 200 basis functions are used. 
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Here the matrix elements /„,'.)„" (£) of the scattering amplitude 
describe the transition from the initial channel with the trans- 
verse energy Ej^'''"^' = /![©! (ni + 5) + a)2(«2 + 5)] and the rel- 



ative longitudinal momentum hk„^ „ 



2pi{E-E 
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y/2pLE\\ (£11 being the relative longitudinal collision energy 
and E the total energy) to the final open channel {n'y^n'2) 

(n' n' ) 

with £ =£'j^'' ^ +EL 0„,,„2(jc,y) are the eigenfunctions of 
the 2D harmonic oscillator corresponding to the eigenvalues 
£]"' "^ rpjjg j^ffgj- ^j-g degenerate with respect to the quantum 
number « = «i +«2 in an isotropic trap (i)\ — oyi — Ox and 



III. RESULTS AND DISCUSSION 
A. Multichannel scattering in isotropic waveguides 

Let us first analyze atomic collisions in a transversely 
isotropic waveguide ©1 = fib = (0^. The partial transmission 
coefficients describing the transmission in the ground Too and 
in the first excited 720 and T{)2 scattering channels are shown 
in Fig[T]for different longitudinal collision energies En. Ob- 
viously the position of the CIR is stable with respect to vari- 
ations of the energy for 10^''^ < E\\/E^ < 10^^. Indeed, this 
position is very close to the value a^/cis = 1 .4603... predicted 
in the zero-energy limit with the pseudo-potential approach 
fill and has been confirmed in subsequent numerical compu- 
tations with different interatomic potentials la [iSl Uql . Here- 
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FIG. 2. (color online) Transition probability PQ2{E,a^/as) from the 
initial ground n = iii + 112 = into the transversally excited states 
)j = «i + 112 = 2 for a few values of a^/as- 



after we define the position of the CIR in the ground state 
scattering channel as the minimum of the transmission co- 
efficient TQo{a±/as,E) which coincides with the zero of the 
imaginary part of the scattering amplitude /oo(flx/fl.s-,£') in 



the region Ej_ 



<E<E 



(2) 



Ol|l[ll[il]. The minimum of the 
transmission coefficient T2o{aj_/as) = TQ2{a±/as) in the first 
excited state (the position of the CIR in the excited scattering 
state) is at a nonzero value and the corresponding transmis- 
sion valley is therefore much less pronounced than in the case 
of 7()o(flx/fl.s) for the ground state. Moreover, the width of 
the CIR differs considerably for the two cases and unlike the 
ground state the position of the CIR in the excited state is 
strongly dependent on the energy Eh. These facts, as we will 
see below, are of crucial importance for the analysis of the 
region near the CIR. Importantly, according to (1) the total 
transmission coefficient is determined in a broad neighbour- 
hood of the resonant region a^ja^ ^ 1.4603... by the partial 
coefficients 72o and T()2 (if W20 and W02 are large enough) 



T{a^/as) w W2o72o(a±/ai) +^027b2(a±/as 



(5) 



which is due to the near-zero values of 7oo(flx/fl.s) around its 
minimum (see FiglT]i. 

A note is in order here. The significant occupation of 
transversally excited states does not necessarily arise due to 
a finite temperature in thermal equilibrium where the occu- 
pation probability is determined by the corresponding Boltz- 
mann distribution. Instead, the preparation of the ultracold 
atomic ensemble in the waveguide itself can lead to an occu- 
pation of excited states. This could e.g. be a nonadiabatic 
loading process into the waveguide. The resulting nonequilib- 
rium state would show a redistribution of its energy between 
the longitudinal and transversal degrees of freedom and might 
eventually thermalize, depending on the number of atoms and 
integrability aspects of the underlying system Ill9ll . 

Atomic collisions near the CIR could possibly also lead to 
an increase of the excited state population. Indeed, the in- 



crease of the atomic loss and heating near a CIR has been at- 
tributed I1IO1I to inelastic three-body collisions, which lead to 
the formation of molecules while transferring the molecular 
binding energy and the following energy release due to possi- 
ble deexcitation of rovibrationally excited molecules to the ki- 
netic energy of the center of mass motion of the molecule and 
in particular to the escaping third atom. Subsequent collisions 
of this third atom can lead to transverse excitations if its ki- 
netic energy exceeds the threshold value 2hoi±^ ~ e\^ — E^ . 
The binding energy of the most weakly bound molecular state 
in the waveguide (see Fig. 2 in lllSll ) exceeds already this 
threshold energy. Using the calculated scattering amplitudes 

jnl'^l we have evaluated the transition probabilities P„„i la] 
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where « = « 1 + «2 and «' 



«2 • The calculated proba- 



bilities P()2(£^,fl±/fls) are presented in Figl2]for a few values 
of flx/fl.v in the region Zs^ < E < E]^ i.e. between the first 
and second excited channel thresholds. Since the probabil- 
ity Po2{E,a±/as) approaches to the values ^ 0.2 — 0.4 one 
might consider collisions with the atoms accelerated due to 
molecules formation near CIR as a possible mechanism for 
the emergence of considerable populations of the first excited 
states n = 2. 

Note also that hereafter we do not address the transmission 
coefficients !;,,„, with odd quantum numbers «i =n2 = 1,3, ... 
because these states are not coupled with the ground and ex- 
cited states possessing even «i and «2. 



B. Multichannel scattering in anisotropic waveguides 

We now analyze the multichannel scattering amplitude and 
the corresponding transmission coefficients for anisotropic 
waveguides for different values of the ratio 01/0)2 7^ 1. In 
Figl3] we present the calculated partial transmission coeffi- 
cients T„^„^{a±/as) as a function of a±/as for fixed values 
of Oi/oh and En. Comparing this result with the isotropic 
case given in Figll]we see that the anisotropy does not change 
the overall behavior of the coefficients Tqo and T20 in the re- 
gion near the CIR. However, the anisotropy splits the excited 
states with « = «i +712 = 2 into two components and changes 
the Tq2 coefficient dramatically by splitting the well of the 
transmission curve. As it is demonstrated in FigH] the ef- 
fect of the splitting of the minimum of the partial transmis- 
sion coefficient To2{a±/as) can also be observed in the total 
transmission coefficient T{a±/as) (1) at 5% of the relative 
population W2/W0 of the first exited states and persists with 
varying 0)1/0)2 and En. Fig |5] demonstrates the dependence 
of the splitting on the relative population W2/W0 calculated in 
the region W2/W0 ;< 0.1 compatible with the experiment llioll . 
The effect of the splitting of the transmission coefficient is 
enhanced with increasing excited state population. 

In the computations, the colliding energy En was chosen in 
agreement with the conditions of the experiment lUOll . where 




FIG. 3. (color online) Partial transmission coefficients 

7^,,„2(a^/ai,£||) in an anisotropic waveguide tt)i/o>2 = 1.05 
as functions of a^ja^ calculated for a near-threshold collision 
energy E\\/e\ = 5 x lO^-'. 



En, even at maximal heating, has remained below Rb x 2>QnK, 

i.e. E\\/Ef^ < 30nK/60QnK = 5 x 10"^ jj should be noted 
that the pronounced shift of the minimum of the T coefficients 
with varying E\\ is due to the considerable dependence of the 
partial coefficients T02 and T2Q on Eh (see FigdJ. This ef- 
fect might be responsible for the shift of the maximum of 
the atom loss in the experiment llioll in the direction of in- 
creasing values for a^. Using this assumption we can fix 
Eu ±2 2.5 x lO^'^E]^ ' as being closest to the experimental con- 
ditions lUOll by choosing from the calculated r(flx/fl.s,£||) 
curves for 0)1/0)2 = 1 the one whose position of the minimum 
coincides more close with the point of maximal atomic loss in 
the experiment. 



C. Diatomic weakly-bound and resonant states in anisotropic 
harmonic waveguides 

To clarify the origin of the splitting of the partial coefficient 
T02 we have calculated the spectrum of the near-threshold 
bound state and resonant states of the atomic dimer in the 
confining trap as a function of a±/as (see Fig|6l). The reso- 
nant energies E,-{a±/as) are determined by the position of the 
minimum of the partial coefficient 7bo(fl±/fl.v,£), where the 
two-body total energy E was varied between the thresholds 
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Oh) and the 
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corresponding to the ground state E^""' = | (o)i 

second excited states Ej'^'"^ = h[(Oi{ni 
with « = «i +n2 = 4. These resonant states were defined in 
Ig] as CIRs with non-zero energies. The binding energy Eg 
of the atomic dimer in the harmonic trap with respect to the 
ground state threshold E^ has been calculated by solving the 
corresponding eigenvalue problem. 

In analyzing the results it is important to note the different 
dependence of the position of the CIR obtained for the differ- 




FIG. 4. (color online) Total transmission coefficients T{a^/as,E«) 
in the region a^/os ~ 1.4603... of the CIR for the isotropic case 
(Oi/cDi = 1 as well as for anisotropic waveguides with coi/cDi = 1.05 
and 1.1 calculated for W2/W0 = 0.05. For 0)1/0)2 7^ 1 the splitting of 
the minimum of the transmission coefficient can be observed. 



ent transmission curves T(io{a±/as),To2 = T2o{a±/as) in FigH] 
on the collisional energy E\\. Since the resonant energy E,- be- 
longing to the CIR in the ground state strongly changes with 
varying a^jus near the value a^ja^ ~ 1.4603... (see Figure|6] 
where a steep descent can be observed with increasing a^ja^ 
the dependence of 7bo(ax/a.v,£) on the energy En near the 
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FIG. 5. (color online) Total transmission coefficients 

T{a^/as,W2/WQ) in the region a^/a^ ~ 1.4603... of the CIR 
for anisotropic case tt)i/ttb = 1.1 as a function of relative population 
W2/lVo of the first excited manifold n = 2. The splitting of the 
minimum of the transmission coefficient is still present even for only 
one percent population of the excited states. 



minimum position is very weak (see Fig[TJ. Opposite to this, 
the dependence of the resonant energy Er for the first excited 
states on a^/cis is much smoother (see uppermost black curve 
with full dots in Fig|6]l. This leads to a considerably stronger 
dependence of the minimum of 72o and Tqi on a^/us on E\\ 
(see Figll). 



In the anisotropic waveguide, the resonant curve E 

,(2) 



(0) 



< 



Er{a^/as) < E^ of the CIR of the ground state splits into 
two components which are shown in Fig|6] The energeti- 
cally lower (0,2)-component qualitatively repeats the behav- 
ior of the resonant energy curve (CIR) of the isotropic case. 
This is why the anisotropy causes only a limited quantita- 
tive change for the Tqq coefficient. However, the behavior 
of the energetically upper (2,0)-component is different (see 
Fig|6]l. The resonance curve E, of this component is flat 
for the complete parameter region of a±/as which leads to 
a strong change with respect to the energy dependence of the 
Tq2 coefficient, including in particular a strong change of the 
corresponding positions of the minima (maxima) with respect 
to a±/as ■ The behavior of the T20 coefficient is determined by 
the energetically lowest resonant curve Er emerging from the 
« = 4 threshold in the isotropic case which is deformed only 
slightly by the anisotropy. This is why the difference between 
72o(wi/fi)2 7^ 1) and Tio (0)1/0)2 = 1) is significantly smaller 
when compared to the case of the T02 coefficient. 

The above model could potentially also explain the appear- 
ance of additional CIRs with further increase of the anisotropy 
COi/oh 7^ 1 as seen in uUA . Increasing 0)1/0)2 7^ 1 the en- 
ergetical distance between the sub-levels characterized by 
n = 2,4,... decreases. This leads to a considerable increase 
of the populations W,^,„2 of higher excited states and in par- 
ticular to an increase of the contribution of these states to the 
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FIG. 6. (color online) Illustration of the spectrum of the atomic 
dimer in harmonic isotropic (Bi /ft)2 = 1 and anisotropic ft)i /ft)2 = 1.2 
waveguides as a function of a^/cis. Below the continuum thresh- 
old E^ the calculated binding energy of the weakly-bound state is 

shown. Between the thresholds E^ and E^ the calculated reso- 
nant energies £,- were determined by the minimum of the transmis- 
sion coefficient Too{a^/as,E). Between the Ej^ and E^ thresh- 
olds the resonant energy Er coincides with the location of the zero of 
/m{/oo(flx/fli)}. Bold solid curves with the full dots correspond to 
the isotropic case a)i /a)2 = 1 and the thin curves with open circles to 
the anisotropic case 0)1/0)2 = 1-2. Indices («i,«2) label the splitted 
sub-levels of the threshold energies E"''"~{a)i/oh = 1.2) of the ex- 
cited states with « = «i +«2- For the second excited threshold with 
n =4 only the lowest sub-level (0,4) is shown. 



total transmission coefficient T (see eq.(l)). 



IV. MECHANISM OF THE CIR SPLITTING IN 
ANISOTROPIC HARMONIC WAVEGUIDE 

The key for the understanding of the mechanism of the 
splitting of the CIR under the action of an anisotropic har- 
monic trap is the diatomic spectrum of the weakly-bound and 
resonant states in the harmonic waveguide given in Fig. 6. It 
also explains why the previous considerations in the zero- 
energy limit near the ground state threshold did not provide 
any splitting of the CIR defined as the singularity of gio = 
limE^^^o{kooRefo/Imfo) MM- 

Actually, so far it was implicitly supposed that the reso- 
nant energy curve £, behaves linearly as it crosses the ground 
state threshold « = for aj_/as = 1 .4603 and becomes subse- 
quently a weakly-bound state for a±/as > 1.4603 (see, for 
example. Fig. 2 in lUSll ). Therefore, it was natural to ex- 
pect that in an anisotropic waveguide 0)i — 0)2 = Ao) -—>■ 
leading to the splitting of the first excited threshold E)^ ' — 

E^' = hAco (see Fig. 6 in the present paper and Fig. 1(b) 
in llOiD the resonant curve Er will also split into two com- 
ponents Er ■ and £,- ' crossing the ground state thresh- 
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old at the points a^/a) ' and a^/al ' with the separation 
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-,(2.0) 



a^/a\. ' — aj_/as '"' proportional to the threshold splitting 
hAco. These expectations were confirmed in the experiment 
by measuring the distance between the maxima of the atomic 
loss in the anisotropic waveguide which was interpreted as 

a±/as ' -fl±/4 ■ (see Fig. 3(c) in ifl^ ). 

However, our extension of the calculation of the resonant 
energy Er to the region a±/as > 1 .4603 has shown a strongly 
nonlinear behavior of the curve Er (see Fig. 6) while shifting 
the point where the resonant curve E,- crosses the ground state 
threshold to the value a±/as — )> +°°- This means that at the 
point fl, -^ +0 we observe a complete rearrangement of the 
spectrum. The first resonant state Er becomes a new weakly- 
bound state once we cross this point and equivalently for the 
higher excited resonant states which convert into each other. 

This behavior, i.e. the rearrangement of the spectrum of 
the two-body system at the point a^ — > +0, remains in the 
anisotropic waveguide (see Fig. 6). This is why the splitting of 
the singularity of the function gio = limEn~^o{kooRefo/Iinfo) 
was not observed near the point oj/as = 1.4603 in the 
anisotropic harmonic waveguide lUU [l2ll but one can observe 
the splitting of the minimum in the effective transmission co- 
efficient T (eq.(l)) due to the quasi-crossing of the resonant 

curves £,■ '' and Er ' leading to the splitting of the mini- 
mum in the T02 coefficient (see Fig. 3) qualitatively equal to the 
width of the quasi-crossing nearox/flv = 1.4603. This width 
is proportional to the first excited threshold splitting hAco in 
the presence of the anisotropy and in good agreement with the 
experimental value for the splitting of the maxima of the atom 
losslEr 



V. CONCLUSIONS 

Our investigation and following analysis demonstrates that 
multichannel scattering in anisotropic harmonic waveguides 



can lead to a splitting of the confinement-induced resonance. 
A necessary ingredient is a population of at least a few percent 
of the transversally excited states which can certainly occur 
via e.g. a nonadiabatic loading process of the atoms into the 
waveguide. 

There are several ways to improve our suggested model 
which would help further clarifying the behaviour of the split- 
ting effect with varying parameters and depending on the ini- 
tial preparation of the atomic ensemble. First, one would 
have to take into account the velocity distribution (distribu- 
tion over the collision energy E\\)in the longitudinal direction 
1201 . Shifts and splittings of the CIR due to anharmonicities 
of the trap 11131 [l4ll and the influence of the inelastic channel 
of molecule formation have also to be determined. 

It is known that the anharmonicity of the trap couples the 
relative and center-of-mass motion of the colliding atoms 
and can lead to additional anharmonicity-induced resonances 
(AIRs) due to the removal of the degeneracy of the center-of- 
mass and relative motions lll3L [l4ll . However, the AIRs are 
supposed to be much narrower than the CIRs and, as a conse- 
quence, more difficult to be detected experimentally because 
of the relatively weak anharmonic coupling with respect to the 
interatomic interaction JMIlSIl- As a conclusion, the observa- 
tion of these AIRs presumably represents a challenging exper- 
imental problem. Explicitly suppressing the anharmonicity of 
the waveguide could lead to a discrimination between the two 
mechanisms causing a splitting of the CIRs. 
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